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Abstract 

A lattice Delaunay poly tope P is called perfect if its De- 
launay sphere is the only ellipsoid circumscribed about P. 
We present a new algorithm for finding perfect Delaunay 
polytopes. Our method overcomes the major shortcomings 
of the previously used method [Du05]. We have imple- 
mented and used our algorithm for finding perfect Delau- 
nay polytopes in dimensions 6, 7, 8. Our findings lead to 
a new conjecture that sheds light on the structure of lattice 
Delaunay tilings. 



1 Introduction 

Let A be an n-dimensional lattice (n > 0) and let P C 
A <g) R = R™ be a poly tope whose vertex set vert P belongs 
to A. We say that Pisa Delaunay polytope for A if P can be 
circumscribed by a closed ball Sp C A <E> R such that 13 p n 
A = vert P. The ball 23 p (or its boundary) is commonly 
referred to as the Delaunay sphere (or empty sphere) for 
P. (Delaunay [Del] himself attributed the concept of empty 
sphere to Voronoi.) Delaunay polytopes for A form a face- 
to-face tiling of A <g> R called the Delaunay tiling for A. 

One can study the geometry of lattices by comparing 
their Delaunay tilings. Such study was initiated by Voronoi 
[VorII08-09]. As A C A ® 1 = K" is deformed into 
Z n C R n by an affine transformation x i— > A(x), the empty 
spheres circumscribed about the Delaunay polytopes of A 
are deformed into empty ellipsoids circumscribed about the 
A-images of these polytopes. All these ellipsoids have iden- 
tical quadratic parts - indeed they are balls in the metric 
d{x,x') = ||.A -1 (a:) - J 4~ 1 (a; / )||. Thus, the study of 
Delaunay tilings for n-lattices is equivalent to the study 
of Delaunay tilings for Z n with respect to different pos- 
itive definite quadratic forms. Let us denote the Delau- 
nay tiling for Z™ with respect to a positive quadratic form 
Q by Del(Z n , Q). The Delaunay property of an ellipsoid 
£(Q,c,P) = {x E R™ | Q[x - c] < R 2 }, circum- 



scribed about a polytope P, means that the quadratic func- 
tion Q[x — c] — R 2 is zero on vert P and strictly positive on 
Z n \vert P. From now on we will be working with A = Z n , 
unless stated otherwise. 

It is natural to extend the notions of Delaunay polytope 
and tiling to positive semidefmite forms. We refer to an R- 
valued function / on a set S as positive if f(x) > for 
any x E S. Let Q be a quadratic form that is positive on 
R™ and such that the rank of the sublattice ker R Q n Z™ 
is equal to the dimension of kerR Q. Then R n is tiled by 
unbounded ?i-dimensional polyhedra, which are Delaunay 
with respect to Q: each polyhedron P from this family is 
circumscribed by an elliptic cylinder £p, whose interior is 
free of lattice points, so that P fl Z™ = £p n Z™. Further- 
more, P is the direct affine product of a Delaunay polytope 
P for a sublattice A C Z" of rank r = dim kerR Q and 
an affine (n — r)-subspace of R™. The degenerate Delau- 
nay "ellipsoid" £p for an unbounded polyhedron P is the 
direct affine product of the Delaunay ellipsoid for P and 
L; we will be using 'ellipsoid' for both bounded and un- 
bounded ellipsoids. For example, R n is tiled by the unit 
slabs Ui = {x | i < xi < i + 1} where i E Z. Each 
unit slab is a Delaunay polyhedron with respect to quadratic 
form x\\ furthermore, each Ui coincides with its Delaunay 
ellipsoid £^. Following a common convention we will be 
using the word 'polytope' only for bounded polyhedra. 

Let A be a lattice of rank n and let P C A® R be a lattice 
polyhedron, i.e., the convex hull of a subset of A. Then P is 
called perfect if there is an n-ellipsoid (possibly degenerate) 
circumscribed about P and this ellipsoid is unique. Perfect 
Delaunay polyhedra are also called extreme (e.g. [Du05]). 
Perfect Delaunay polytopes are rare in small dimensions, 
e.g., for n < 6 there are only three such polytopes - for 
n = 0, [0, 1] for n = 1, and Gosset's semiregular polytope 
22i for n — 6. The previous method [Du05] for finding per- 
fect Delaunay polytopes was based on an unproven conjec- 
ture that every perfect Delaunay polytope is basic. A lattice 
polytope P is called basic if there exist vq, . . . , v n G vert P 
such that every v E vert P can be written as v = X)T=o ^ iVi 



where Ai, . . . , A„ + i G Z and 1 = J2?=o ^ We know that 
there exist non-basic Delaunay polytopes in higher dimen- 
sions (see [DG07]) and we cannot rule out the existence of 
non-basic perfect Delaunay polytopes. 

The perfection property of a Delaunay polytope P with 
ellipsoid £(Q, c, p) = {x e R" | Q[x - c] < p 2 } amounts 
to that any quadratic function that vanishes on vert P is of 
the form a(Q[x — c] — p 2 ) where a 6 R. A real-valued 
quadratic function F on R™ is called perfect if armin F = 
min{F(z) | z e Z™} > and conv{z e Z™ | F(z) = 
armin F} is a perfect Delaunay polyhedron. The ellipsoid 
{x £ R™ I P(:e) < arming 1 } is also called perfect. 

Erdahl proved that the vertex set of any perfect Delaunay 
polyhedron splits uniquely (up to arithmetic equivalence) 
into the direct affine sum of the vertex set of a perfect Delau- 
nay polytope and a sublattice which is parallel to the kernel 
ofQ[Er92]. 

Theorem 1 [Er92] A polyhedron P 6 Del(Z n , Q) is per- 
fect if and only if 

PC\7L n = {v + z\v e vert D, z E F}, 

where D is a perfect polytope from Del(Z n n affD, Q) 
and T is a submodule of IT 1 such that Z™ is the direct sum 
of modules (Z™ n aff D) - (Z™ n aff D) and T. If{D', V) is 
another pair with these properties, then V = Y and D 1 = 
A(D), where A is an affine automorphism of IT 1 . 

The fundamental importance of perfect Delaunay poly- 
hedra is explained by the following theorem of Erdahl 
[E00]. 

Theorem 2 Let D be a Delaunay polytope for Z". Then 

k 

D = f]P, 

i=l 

where each Pi is a perfect Delaunay polyhedron for Z" with 
respect to some positive form Qj such that dim kern Qj = 
rank(ker R Q 4 n Z n ). 

For S C Z" we define qrank S, the quadratic rank of S, 
as the dimension of the space of quadratic functions on R™ 
that vanish on S. Let P and Pi be two perfect Delaunay 
polytopes for Z™ such that qrank(vert P n vert Pi) = 2. In 
this case we call P and Pi adjacent. If two perfect Delau- 
nay n-polytopes P and P' can be connected by a sequence 
of perfect Delaunay n-polytopes in which every two con- 
secutive members are adjacent, then we will say that P and 
P' belong to the same adjacency component. We developed 
a method for finding an adjacency component for a perfect 
Delaunay polytope. We have found that for each n < 8 
all known perfect Delaunay n-polytopes belong to the same 
adjacency component. This finding makes compelling the 
following conjecture. 



Conjecture 1 For any n G N all perfect Delaunay n- 
polytopes belong to the same adjacency component. 

2 Space of quadratic functions 

Let us denote by Sym(n) the space of real symmet- 
ric n x n matrices; by interpreting an element of Sym(n) 
as a Gram matrix, we can regard Sym(n) as the space of 
quadratic forms with real coefficients. Denote by Q(n) the 
linear space of quadratic functions on W 1 and by Q (n) C 
Q(n) the subspace of functions with zero constant term. 
Since a quadratic function can be represented uniquely as 
the sum of a quadratic form, a linear functional, and a con- 
stant, it is convenient to introduce the projection operators 

Quad : Q(n) Sym(n), Lin : Q(n) -> M™*, 

Const : Q(n) R. 

For two quadratic forms with Gram matrices A and B we 
define the dot product on Sym(n) as trace(AB). For linear 
functions defined by covectors a and b the dot product is just 
a ■ b. The dot product on Qo(n) is defined as the direct sum 
of the dot products on Sym(n) and M ra *. 

The main idea of this paper is in interpreting quadratic 
functions on R" as elements of Q (n)*, the dual of Qo(n). 
There is a natural correspondence between ellipsoids in M™ 
and closed (affine) halfspaces of Qo(n). Namely, if £ = 
{x e M™ | Q[x - c] < p 2 }, then the corresponding halfs- 
pace H E is {X e Q (n) | X ■ F > p 2 - Q[c]}, where F is 
the quadratic function defined by F(x) = Q[x] — 2Q(x, c). 

Let D be a map from M™ into Qo(n) defined in the matrix 
notation by 

!D '. u i ► (x i > x T (uu T )x + u T x), 

where x and u are treated as column vectors. Obviously, 
D takes an integer vector to a quadratic function with in- 
teger Gram matrix and integer linear part; such quadratic 
function is called classically integer. The map D resembles 
the Voronoi map V : R™ — ► Sym(n) that takes a vector u 
to the quadratic form with Gram matrix uu T (see [RB79] 
for details). Thus, we have D(u) = V(u) + u*, where u* 
is linear functional dual to it. The map V can be seen as 
the quadratic Veronese map from R" to Sym(n), although 
in contemporary literature the Veronese map is usually de- 
fined in the projective setup. We call an ellipsoid in R™ 
empty if its interior is free of points of Z". If £ = {x £ 
R" | Q[x-c] < p 2 } is an empty ellipsoid, then Hg contains 
allofD(Z"). Thus, T> (Z") c f] H E . The right 

{£ is empty} 

hand side is the intersection of infinitely many halfspaces, 
which can be replaced by the intersection of only those half- 
spaces whose boundaries are completely determined by the 
elements of D(Z") that lie on them. Throughout the paper 



we use conv S to denote the convex hull of a set S C W 1 
and aff S to denote the minimal affine subspace containing 
S. 

We define the Erdahl polyhedron E(n) as the intersec- 
tion of closed halfspaces Hg such that £ is empty and 

rank(<9iie n D (Z" ) ) = dim Q (n) . 

Theorem 3 TTie set D(Z n ) coincides with X>(R n ) n<9E(n). 
Furthermore, 

conv'D(Z n ) = E(n). 

Proof. Each z G Z™ belongs to the boundary of an empty 
degenerate ellipsoid 

£ = {a; G K" | z • ei < x • z < z • ei + 1}. 

It is easy to check that <9£ is the only quadratic surface pass- 
ing through Z™ (~l <9£ (or see [Er92] for a proof). Thus, 
rankaHg n D(Z") = dimQ (n) and D(Z") C <9E(n). 

Let X G D(K n ) fl 9E(n). Then there is x G R™ with 
X = T)(x). If x i Z" (for x G Z" see above), then 
x = u + x', where it G Z" and x' £ [0, 1]". Without loss 
of generality assume that x[ £ Z. Then x 6 int £, where 
£ = {x | u\ < x\ < ui + 1}, which means X ^ iJg, 
contradicting our choice of X. Thus, the only points of R™ 
that are mapped by CD on 9 E(n) are elements of Z" and the 
first claim of the theorem is proven. 

Suppose X G E(n) and let us show X G convD(Z"). 
It is enough to prove this implication for X G 9E(n). If 
X 6 9E(n), then X lies on some Hg, where £ is empty 
and completely determined by elements of Z" that lie on its 
boundary. If X ^ conv CD(£ n Z n ), then there is a facet C 
of conv D(£ n Z") such that X and the relative interior of 
convD(£ n Z n ) lie in the different halfspaces of Hg with 
respect to aff C. Let /e(x) > be an affine inequality 
defining Hg and let gc{x) — be an affine equation of the 
hyperplane passing through the point I (where I : x ^ 
x T x) and aff C such that gc{X) < 0. The equation of any 
hyperplane in Qo(n) passing through aff C can be written 
as /g (x) + 9gc (x) = for some pel. Let us define 

6 m - sup{0 G R | Vz G Z" / e (D(z))+^c(2)(z)) > 0} 

We claim that m > and there exists u e Z n \ E such 
that 

/£(D(u))+0 rofl c(D(u)) =0. 

The proof of these claims (which we omit due to the space 
limitations) is based on standard techniques of geometry of 
numbers and follows the line of argument used by Voronoi 
in his first memoir [VorCol, Pages 177-179]. □ 

If an ellipsoid £ is empty and dH^ contains 
dimQo(n) + 1 affinely independent points, then £ is 
uniquely determined by the points of Z™ that lie on its 



boundary: in this case £ is called a perfect ellipsoid for 
lattice Z™. Perfect ellipsoids were introduced by Erdahl 
[Er75,Er92]. Thus, 

E(n) = conv'D(Z™) = f] H e . 

{E is perfect} 

Note that E(n) is not a polyhedron in the sense of lin- 
ear programming, where the number of constrains is al- 
ways assumed to be finite. We will refer to the faces 
of E(n) of dimension dimQo(n) — 1 as facets and the 
facets of dimension dimQ (n) — 1 as faces. The facets 
of E(n) correspond to perfect Delaunay polyhedra. The 
bounded facets of E(n) correspond to perfect Delaunay 
poly topes. Two perfect Delaunay poly topes are adja- 
cent if the corresponding facets of E(n) share a bounded 
ridge. Faces of E(n) correspond to Delaunay polyhedra 
- bounded faces to bounded Delaunay polyhedra (poly- 
topes) and unbounded faces to unbounded polyhedra. There 
is a great deal of analogy between E(n) and Voronoi's 
polyhedron II(n), introduced by Venkov [Ven40] (see also 
[RB79]). Recall that II(n) is defined as the convex hull 
of {V(p) | p G Z n and g.c.d.(pi, . . . ,p n ) = 1}, where 
V : R™ -> Sym(n,R) is the Voronoi map. The facets of 
II(n) are defined by closed halfspaces corresponding to per- 
fect forms, which were studied by Voronoi [VorI08] (part I). 
The bounded facets of IT(n) correspond to positive definite 
perfect forms. 

2.1 Geometry of E(n) 

Denote by Aff n (Z) the group of affine automorphisms 
of Z™, i.e. the group of transformations of the form A(z) = 
Lz + t, where L G GL n (Z) and t G Z™. The action of 
Aff n (Z) on R" can be naturally lifted to Q (n) by 

F i — ► {^x i — ► F(A~ 1 (x))}. 

The group A//„(Z) acts on E(n) in a way somewhat 
similar to that of GL n (Z) acting on II(n). Subsets V and 
V' of R™ are called arithmetically equivalent if there exists 
A G Aff n (Z) such that A(V) = V. Obviously, arith- 
metic equivalence preserves properties of ellipsoids such as 
the Delaunay property, emptiness, and perfection. Since 
there are only finitely many arithmetically distinct Delau- 
nay polytopes in each dimension (e.g. [DL97]), the bound- 
ary of E(n) has finitely many distinct arithmetic types of 
faces. In fact, the definition of perfect ellipsoid implies that 
arithmetically equivalent perfect ellipsoids are isometric. 

There are beautiful connections between the polytope 
E(n) and Delaunay tilings of Z n . The projection Lin : 
Qo(n) — > R™ maps the vertices of 3E(n) onto the points 
of Z n . The projection of each face of E(n) is a Delaunay 
polyhedron in Del{Z n , Q) for some positive quadratic form 
Q. In particular, the projections of facets of E(n) are perfect 
Delaunay polyhedra. 



3 Algorithm 

In this paper we present a practical algorithm that finds 
all perfect Delaunay polytopes that belong to the adja- 
cency component of a known n-dimensional perfect De- 
launay polytope. Our algorithm is best explained geomet- 
rically in terms of the geometry of Qo(n), although it is 
easier to implement it in terms of Q(n) by representing the 
closed halfspace iJg corresponding to an empty ellipsoid 
£ = {x e R" | F(x) < 0} by the ray R+F in Q(n). 
In this dual interpretation we consider the convex hull C(n) 
in Q(n) of all rays corresponding to all empty ellipsoids of 
E(n). Each extreme ray of the cone Q(n) is of the form 
RP, where P is a perfect quadratic function. Thus, the ad- 
jacency between the facets of E(n) corresponds to the adja- 
cency between the extreme rays of the cone 6(n), i.e., facets 
of E(n) determined by perfect functions F and F' are ad- 
jacent if an only if the rays M.F and M.F' share a common 
2-face of the cone C(n). The convenience of this represen- 
tation for computing is much due to its homogeneity, that is, 
affine transformations of R™ induce linear transformations 
on Q(n). 

We record our knowledge of the adjacency compo- 
nent under investigation in the adjacency graph G(V, E), 
where V is the set of arithmetic types of perfect Delau- 
nay polytopes and £ is the set of arithmetic types of pairs 
(P, P'), where P and P' are perfect Delaunay polytopes 
and q ran k(vert P n vert P') = 2. Equivalently, one can 
think of V as of a set of inequivalent facets of E(n) and £ 
as the a set of inequivalent ridges of E(n). In the following 
subsection we describe the basic step of the algorithm. 

3.1 Step of the Algorithm 

Let P £ DeZ(Z™, Q) be a perfect Delaunay polytope and 
let £ = {xel" F(x) <a,Fe Q (n)} be its empty 
circumscribed ellipsoid. We regard P as a vertex of the 
adjacency graph G(V, E). At each moment the algorithm 
is looking at a particular vertex of this graph. First we find 
a subset S of vertP with qrankS" = 2. Let Hg = {X G 
Qo(n) | F ■ X > a} and let G ■ X = (3 be the equation of 
the hyperplane in Q (n) passing through the point I (where 
I{x) = x T x) and aff 2)(S) such that G ■ D{v) < (3 for 
all v £ vert P. The equation of any hyperplane passing 
through aff D (S) can be written as F ■ X + pG ■ X = a + pf3 
for some pgR. Let 

p m = sup{p e R | Vz e Z" F-D(z)+pG-T>(z) > a+p(3}. 

In situations like this it is often said that the hyperplanes 
{K{p)}, where 

9t(p) = {X\F-X + pG-X = a + p(3}, 



hinge on the ridge aff D(S) n E(n) of the surface 9E(n) 
and that p is the hinge parameter. 

It can be shown (see Theorem 3) there exists u e Z n \P 
such that 

F ■ -D(u) + Pm G ■ D(u) >a + Pm f3. 

The search for msZ"\P such that F ■ D(u) + p m G ■ 
CD (it) > a + p m f3 can be interpreted as continuous rotation 
of the hyperplane CK(p) from the initial position at p = 
to the final position at p = p m (see Figure 1). For small 
values of p the hyperplane 'K(p) intersects with E(n) only 
over the ridge aff D{S) fl E(n). When p reaches the value 
of p m the rotational motion of the hyperplane is stopped by 
the point CD (u). S and u define a perfect quadratic function 
and corresponding perfect Delaunay polyhedron. If the new 
polyhedron is a polytope, we check whether it is arithmeti- 
cally equivalent to any of the already discovered polytopes. 
If it is a polytope distinct from the previously discovered 
ones, we add it to the list of perfect n-polytopes and up- 
date the adjacency graph. This procedure is similar to the 
one used by Voronoi in the determination of perfect forms 
in small dimensions. He referred to this procedure as the 
method of continuous variation of parameters. The geomet- 
ric interpretation of Voronoi's method as that of hinging hy- 
perplanes was given by Venkov [Ven40]. Later this method 
was rediscovered in the context of polytopes by [CK70] and 
dubbed as "gift wrapping method". 

The procedures described above paragraph are repeated 
for each arithmetic class of subsets of vert P of quadratic 
rank 2. When all such subsets are exhausted, we move to 
another vertex of the adjacency graph G(V, E) 

3.1.1 Finding p m and u 

Let S C vert P and let qrank S = 2. Using some heuristic 
we pick some z £ vert P with G ■ D(z) > (i and construct 
an ellipsoid £ through S and z. We find its center c and 
then look for the closest lattice point to c in the metric de- 
fined by £. This test can be done efficiently for n < 9 using 
the program Lattice-CVP by Dutour (see [LCVP]). If 
the closest lattice points happens to be at the same distance 
from cas z and S, then we declare Z™ n £ the vertex set 
of a perfect Delaunay polyhedron. If the interior of £ con- 
tains a lattice point z', then we abandon z and repeat the 
computation for S and z', etc. 

3.2 Using Symmetries in computation 

Our algorithm would be impractical if we failed to use 
symmetries in an efficient way. Two isomorphism problems 
had to be addressed. The first is the problem of checking 
whether two perfect Delaunay polytopes are arithmetically 
equivalent. The second problem is finding all arithmetically 



Figure 1. Going from facet conv{(l, 1), (4,2)} to facet conv{(4,2), (9,3)}. 



inequivalent subsets S of vertP with qrank5 = 2. Algo- 
rithms for these problems have been implemented in GAP 
(some of them are available in [DuPol]) and rely on the use 
of the program nauty [McKay]. See also [Du05]. 

3.3 Results 

The graph G(V, E) constructed by the algorithm en- 
codes the adjacency pattern for bounded facets of E(n). 
More generally, denote by G(V, E) the graph whose ver- 
tices are the arithmetic types of facets of E(n) and whose 
edges are arithmetic types of pairs of facets sharing a com- 
mon ridge. As of now this graph is completely known only 
for n < 6. For n = 6 G(V, E) has two vertices, which 
correspond to the Gosset 6-polytope 22i and the unit slab 
U. The quadratic form for the Gosset 6-polytope is E 6 and 
that for the unit slab is a rank one form (see Figure 2). 

For n = 7 the discovered adjacency component of 
G(V, E) has 4 vertices and that of G(V, E) has two ver- 
tices (see Figure 3). The latter two vertices correspond to 
the Gosset 7-polytope 32i and a polytope with 35 vertices 
discovered earlier by Erdahl and Rybnikov (see [ErRyb]). 

For n = 8 we have determined an adjacency component 
of the restricted graph G(V, E). Below is the adjacency list 




u 

Figure 2. n = 6. left: G(V,E), right: G(V, E). 



of the conjectured G(V, E) in GAP format. Note that the 

graph has loops and multiple edges. 

1:[1,2] 

2: [2, 16], [2, 27], [2, 8], [2, 10], [2, 22], [2, 4], [2, 5], [2, 13], 
[2, 7], [2, 6], [2, 3], [2, 14], [2, 12], [2, 19], [2, 9], [2, 18], [2, 8], 
[2, 11], [2, 15], [2, 6], [2, 11], [2, 2], [2, 17], [2, 1] 
3: [3, 20], [3, 13], [3, 12], [3, 3], [3, 11], [3, 4], [3, 14], [3, 5], 
[3, 15], [3, 10], [3, 6], [3, 2] 

4: [4, 5], [4, 6], [4, 10], [4, 22], [4, 3], [4, 2], [4, 8], [4, 14], [4, 20], 
[4, 19] 

5: [5, 9], [5, 21], [5, 6], [5, 10], [5, 5], [5, 9], [5, 6], [5, 22], [5, 8], 
[5, 20], [5, 4], [5, 3], [5, 2] 

6: [6, 22], [6, 10], [6, 5], [6, 4], [6, 10], [6, 3], [6, 2], [6, 9], [6, 6], 
[6, 5], [6, 21], [6, 8], [6, 12], [6, 15], [6, 24], [6, 6], [6, 7], [6, 11], 
[6, 2], [6, 15] 

7: [7, 7], [7, 12], [7, 21], [7, 22], [7, 9], [7, 24], [7, 19], [7, 8], 




35-tope 



Figure 3. n = 7. left: conjectured G(V, E), right: conjectured G(V, E). 



[7, 6], [7, 10], [7, 2] 

8: [8, 22], [8, 2], [8, 27], [8, 8], [8, 16], [8, 8], [8, 10], [8, 20], 
[8, 9], [8, 6], [8, 2], [8, 5], [8, 8], [8, 4], [8, 13], [8, 12], [8, 7] 
9: [9, 8], [9, 6], [9, 5], [9, 22], [9, 5], [9, 7], [9, 2], [9, 10], [9, 19], 
[9,21], [9, 23], [9,20], [9, 22] 

10: [10, 6], [10, 5], [10, 15], [10, 24], [10, 10], [10, 22], [10, 9], 
[10, 12], [10, 21], [10, 7], [10, 10], [10, 10], [10, 3], [10, 2], 
[10, 20], [10, 8], [10, 4], [10, 26], [10, 6], [10, 12], [10, 13], 
[10, 16], [10, 16] 

11: [11,6], [11,3], [11,2], [11,2] 

12: [12, 10], [12, 20], [12, 3], [12, 8], [12, 10], [12, 2], [12, 7], 
[12, 6], [12, 21] 

13: [13, 10], [13, 20], [13, 3], [13, 8], [13, 2] 

14: [14, 4], [14, 18], [14, 3], [14, 2] 

15: [15, 2], [15, 6], [15, 6], [15, 3], [15, 10], [15, 21] 

16: [16, 8], [16, 2], [16, 27], [16, 20], [16, 10], [16, 10], 

17: [17, 2] 

18: [18, 19], [18, 2], [18, 14] 

19: [19, 9], [19, 7], [19, 2], [19, 4], [19, 20], [19, 18], [19, 25] 
20: [20, 16], [20, 22], [20, 9], [20, 5], [20, 19], [20, 4], [20, 20], 
[20, 12], [20, 3], [20, 13], [20, 10], [20, 8] 

21: [21, 7], [21, 12], [21, 24], [21, 21], [21, 22], [21, 15], [21, 5], 
[21, 9], [21, 23], [21, 26], [21, 6], [21, 10] 

22: [22, 22], [22, 2], [22, 8], [22, 6], [22, 27], [22, 9], [22, 7], 

[22, 4], [22, 5], [22, 10], [22, 9], [22, 22], [22, 21], [22, 20], 

23: [23, 9], [23, 21], [23, 25] 

24: [24, 6], [24, 21], [24, 7], [24, 10] 

25: [25, 23], [25, 19] 

26: [26, 10], [26, 21] 

27: [27, 8], [27, 22], [27, 16], [27, 2] 

The numbers of vertices correspond to the numbers of 
polytopes in [DuErRy] where a complete analysis of exist- 
ing data on perfect Delaunay polyhedra for n < 8 is given. 

In practice the algorithm often encounters perfect ellip- 
soids equivalent to the unit slab, i.e., to the set < x\ < 1. 
For n G {6, 7, 8} we know that any perfect polyhedron ad- 
jacent to the unit slab is either a unit slab or the product of 
Gosset's 6-polytope 2 2 i and M n ~ 6 . However, for n > 6 
other unbounded perfect polyhedra appear, such as e.g. the 
product of 22i and M 2 for n = 8. We are not able, at 
the present, to comprehensively handle all subsets S with 
q rankS 1 = 2 for such polyhedra. That is why we cannot 
formally claim that our results for n = 7, 8 are complete. 



4 Discussion 

The new method has many advantages over the one of 
[Du05]: 

1. Unlike the previous methods (see e.g. [Du05]), the 
new method uses the full symmetry group of P. In 
particular, the use of the full symmetry group of P 
allows us to use the Recursive Adjacency Decompo- 
sition Method of [BDS07] to terminate computations. 
The termination problem is an important one. Previ- 
ous methods did not have a satisfactory solutions to 
the termination problem. 

2. Previous methods had to select an affine basis for each 
perfect Delaunay polytope. We do not know if it is 
possible to find an affine basis for every perfect De- 
launay polytope. The method of this paper does not 
require this assumption. 

3. Our method is no longer reduced to basic Delaunay 
polytopes. We know that there exist non-basic Delau- 
nay polytopes (see [DG07]) and we cannot exclude the 
possibility that there exist non-basic perfect Delaunay 
polytopes. 

4. The new method has found all presently known 8- 
dimensional perfect Delaunay polytopes. The method 
of [Du05] run in dimension 8 does not find some 
of these polytopes - they were found as sections of 
higher-dimensional perfect Delaunay polytopes ob- 
tained earlier by the old method; however, these sec- 
tions were found by a heuristic approach without any 
guarantee of completeness. On the other hand, if Con- 
jecture 1 is true, then the new algorithm has provably 
found all perfect Delaunay n-polytopes for n < 8. We 
expect that running the new algorithm for n = 9,10 
will uncover previously unknown perfect polytopes in 
these dimensions. 

Our method cannot deal at present with unbounded perfect 
Delaunay polyhedra. When a perfect polyhedron P is un- 
bounded, but not equivalent to the unit slab, it is difficult 



to find all equivalence classes of its vertex subsets corre- 
sponding to the ridges of E(n). For this reason we cannot 
guarantee that our algorithm has found all perfect ellipsoids 
in dimensions 7 and 8. 
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